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In a field-theoretical context, we consider the Euclidean 
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the spatial dimensions. We are able to determine the dependence of the transition temperature 
(T c ) for a system described by this model, confined between two parallel planes, as a function of 
the distance (L) separating them. We show that T c is a concave function of L" 1 . We determine a 
minimal separation below which the transition is suppressed. 



PACS number(s): 03.70.+k, ll.10.-z 

In the last few decades, a large amount of work has 
been done on field theoretical models applied to the study 
of critical phenomena. In particular, second-order phase 
transitions have been extensively studied in view of the 
investigation of several material systems. An account on 
the state of the subject and related topics can be found, 
for instance, in Refs. Questions concerning the 

existence of phase transitions may also be raised if one 
considers the behaviour of field theories as a function of 
spatial boundaries. The existence of phase transitions 
would be in this case associated to some spatial param- 
eters describing the breaking of translational invariance, 
for instance, the distance L between planes confining the 
system. Studies of this type have been recently per- 
formed 0, 0| , concerning with the spontaneous sym- 
metry breaking in the A0 4 theory. In particular, if one 
considers the Ginzburg-Landau model confined between 
two parallel planes, which is assumed to describe a film 
of some material, the question of how the critical tem- 
perature depends on the thickness L of the film can be 
raised. 

Studies on confined field theory have been done in the 
literature for a long time. In particular, an analysis of 
the renormalization group in finite size geometries can 
be found in 0, . These studies have been performed 
to take into account boundary effects on scaling laws. In 
another related topic of investigation, there are systems 
that present domain walls as defects, created for instance 
in the process of crystal growth by some prepared cir- 
cumstances. At the level of effective field theories, in 
many cases, this can be modeled by considering a Dirac 
fermionic field whose mass changes sign as it crosses the 
defect, meaning that the domain wall plays the role of 
a critical boundary separating two different states of the 
system EH EH. 

Under the assumption that information about gen- 
eral features of the behaviour of systems undergoing 
phase transitions in absence of external influences (like 
magnetic fields) can be obtained in the approximation 
which neglects gauge field contributions in the Ginzburg- 
Landau model, investigations have been done with an 



approach different from the renormalization group anal- 
ysis. It has been considered the system confined between 
two parallel planes and using the formalism developed in 
Refs. [Til ITq it has been investigated how the critical 
temperature is affected by the presence of boundaries. 
In particular a study has been done on how the critical 
temperature (T V) o f a superconducting film depends on 
its thickness L 0, 0] . In this paper we perform a fur- 
ther step, by considering in the same context an extended 
model, which besides the quartic field self-interaction, a 
sextic one is also present. It is well known that those in- 
teractions, taken together, lead to renormalizable quan- 
tum field theories in three dimensions and that they are 
supposed to describe first-order phase transitions. 

/From our point of view, as in previous publications, 
the system to be studied is a slab of a material of thick- 
ness L, the behaviour of which in the critical region is 
to be derived from a quantum field theory calculation 
of the dependence of the renormalized mass parameter 
on L. We start from the effective potential for the the- 
ory, which is related to the renormalized mass through a 
renormalization condition. This condition, however, re- 
duces considerably the number of relevant Feynman di- 
agrams contributing to the mass renormalization, if one 
wishes to be restricted to first-order terms in both cou- 
pling constants of the model. In fact, just two diagrams 
need to be considered in this approximation: a tadpole 
graph with the 4 coupling (I loop) and a "shoestring" 
graph with the cj) 6 coupling (2 loops)(see Fig. I). No di- 
agram with both couplings occur. The L-dependence 
appears from the treatment of the loop integrals, as the 
material is confined between two plane sheets a distance 
L apart from one another. We therefore take the space di- 
mension orthogonal to the planes as finite, the other two 
being otherwise infinite. This dimension of finite extent 
is treated in the momentum space using the formalism of 
Ref. HI. 

We start by stating the Ginzburg-Landau Hamiltonian 
density in a Euclidean .D-dimensional space, now includ- 
ing both 4 and (j> 6 interactions, in the absence of external 
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fields, given by (in natural units, fi 
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where A and rj are the (renormalized) quartic and sex- 
tic self-coupling constants, with the bare mass given by 
tuq = a(T—Ta), Tq being the bulk transition temperature 
of the material and a > 0. We consider the system con- 
fined between two parallel planes, normal to the x-axis, 
a distance L apart from one another and use Cartesian 
coordinates r = (x, z), where z is a (D — l)-dimensional 
vector, with corresponding momenta k = (fc-^q^q be- 
ing a (D — l)-dimensional vector in momenta space. The 
generating functional of Schwinger functions is written in 
the form 



dx / d D - l zH{\if\,\W(p\) 
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with the field <p(x, z) satisfying the condition of confine- 
ment along the x-axis, tp(x < 0, z) = (p(x > L, z) =const. 
Then the field should have a mixed series-integral Fourier 
representation of the form 
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where w n = 27m/ L and the coefficients c n and 6(q) cor- 
respond respectively to the Fourier series representation 
over x and to the Fourier integral representation over 
the (D — l)-dimensional z-space. The above conditions 
of confinement of the x-dependence of the field to a seg- 
ment of length L allow us to proceed, with respect to 
the x-coordinate, in a manner analogous as is done in 
the imaginary-time Matsubara formalism in field theory 
and, accordingly, the Feynman rules should be modified 
following the prescription 
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We emphasize, however, that we are considering an Eu- 
clidean field theory in D purely spatial dimensions, so we 
are not working in the framework of finite-temperature 
field theory. Here, the temperature is introduced in the 
mass term of the Hamiltonian by means of the usual 
Ginzburg-Landau prescription. 

To continue, we use some one-loop results described in 
[T3. ITU adapted to our present situation. These re- 
sults have been obtained by the concurrent use of dimen- 
sional and zeta- function analytic regularizations, to eval- 
uate formally the integral over the continuous momenta 
and the summation over the frequencies uj n . We get sums 
of polar (L-independent) terms plus L-dependent ana- 
lytic corrections. Renormalized quantities are obtained 
by subtraction of the divergent (polar) terms appearing 
in the quantities obtained by application of the modified 



Feynman rules and dimensional regularization formulas. 
These polar terms are proportional to T-functions hav- 
ing the dimension D in the argument and correspond to 
the introduction of counterterms in the original Hamil- 
tonian density. In order to have a coherent procedure in 
any dimension, those subtractions should be performed 
even for those values of the dimension D for which no 
poles are present. In these cases a finite renormalization 
is performed. 

In principle, the effective potential for systems with 
spontaneous symmetry breaking is obtained, following 
the Coleman- Weinberg analysis |2(j , as an expansion in 
the number of loops in Feynman diagrams. Accordingly, 
to the free propagator and to the no-loop (tree) diagrams 
for both couplings, radiative corrections are added, with 
increasing number of loops. Thus, at the 1-loop approxi- 
mation, we get the infinite series of 1-loop diagrams with 
all numbers of insertions of the <f> vertex (two external 
legs in each vertex), plus the infinite series of 1-loop di- 
agrams with all numbers of insertions of the 4> 6 vertex 
(four external legs in each vertex), plus the infinite series 
of 1-loop diagrams with all kinds of mixed numbers of 
insertions of <fi 4 and <fi e vertices. Analogously, we should 
include all those types of insertions in diagrams with 2 
loops, etc. However, instead of undertaking such a daunt- 
ing computation, even if we restrict ourselves to the low- 
est terms in the loop expansion, we remember that the 
gap equation we are seeking is given by the renormaliza- 
tion condition in which the renormalized squared mass is 
defined as the second derivative of the effective potential 
U(ipo) with respect to the classical field (po, taken at zero 
field, 
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For our purposes, we do not need to consider the renor- 
malization conditions for the interaction coupling con- 
stants, i.e., they may be considered as already renormal- 
ized when they are written in the Hamiltonian above. 
At the 1-loop approximation, the contribution of loops 
with only <fi A vertices to the effective potential is obtained 
directly from 15], as an adaptation of the Coleman- 
Weinberg expression after compactification in one dimen- 
sion, 
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In the above formula, in order to deal with dimension- 
less quantities in the regularization procedure, we have 
introduced parameters c 2 = m 2 /A / n 2 fi 2 , a = (L/i)~ 2 , 
gi = \/8ir 2 and (f>a — (fo/js, where ipo is the normal- 
ized vacuum expectation value of the field (the classical 
field) and fi is a mass scale. The parameter s counts the 
number of vertices on the loop. 
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FIG. 1: Diagrams contributing to the renormalized mass at 
lowest order in the coupling constants. 



It is easily seen that only the s — 1 term contributes to 
the renormalization condition JSJ). It corresponds to the 
tadpole diagram. It is then also clear that all 6 -vertex 
and mixed 4 - and (/> 6 -vertex insertions on the 1-loop di- 
agrams do not contribute when one computes the second 
derivative of similar expressions with respect to the field 
at zero field: only diagrams with two external legs should 
survive. This is impossible for a </> 6 -vertex insertion at 
the 1-loop approximation, therefore the first contribu- 
tion from the <fi 6 coupling must come from a higher-order 
term in the loop expansion. Two-loop diagrams with two 
external legs and only (/> 4 vertices are of second order in 
its coupling constant, and we neglect them, as well as all 
possible diagrams with vertices of mixed type. However, 
the 2-loop shoestring diagram, with only one cj) e vertex 
and two external legs is a first-order (in 77) contribution to 
the effective potential, according to our renormalization 
criterion. 

Therefore the renormalized mass is defined at first- 
order in both coupling constants, by the contributions 
of radiative corrections from only two diagrams: the tad- 
pole and the shoestring diagrams. The tadpole contribu- 
tion reads (putting s = 1 in eq. ©), 

U 1 (<p ,L)=p D ^9i<t>t £ / k2 „„ 2 „ 2 - (7) 



The integral on the D — 1 non-compactified momentum 
variables is performed using the dimensional regulariza- 
tion formula 
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The sum in the above expression may be recognized as 
one of the Epstein-Hurwitz zeta- functions, Zf (^^;a), 
which may be analytically continued to [21| 
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where the K v are Bessel functions of the third kind. The 
tadpole part of the effective potential is then 
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Notice that since we are using dimensional regularization 
techniques, there is implicit in the above formulas a factor 
p^ D in the definition of the coupling constant A. In 
what follows we make explicit this factor, the symbol A 
standing for the dimensionless coupling parameter (which 
coincides with the physical coupling constant in D = 4). 

We now turn to the 2-loop shoestring diagram contri- 
bution to the effective potential, using again the Feyn- 
man rule prescription for the compactified dimension. It 
reads 



U 2 {4>,L)=f D - 2 a^g 2 ct>l 
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or, after subtraction of the polar term coming from the 
first term of Ea.lfTUfl. 
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The renormalized mass with both contributions then 
satisfies an L-dependent generalized Dyson-Schwinger 
equation, 
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where we have introduced the dimensionful coupling con- 
stants (3 = — \p 4 ~ D and p — ■np 6 ^ 20 . 

A first-order transition occurs when all the three min- 
ima of the potential 



V(y Q )=m(LYvl-^t + fyl 
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where m(L) is the renormalized mass defined above, are 
simultaneously on the line V(ipo) = 0. This gives the 
condition 
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FIG. 2: Relative critical temperature v = T c /Tb as a function 
of the inverse thickness x — 1/L, in natural units {m~ l ) with 
To = 4000, a = 1, 13 = 0.1 p= 0.3/4. 

Notice that the value to = is excluded in the above 
condition, for it corresponds to a second-order transition. 
For D = 3, which is the physically interesting situation 
of the system confined between two parallel planes em- 
bedded in a 3-dimensional Euclidean space, the Bessel 
functions entering in the above equations have an ex- 
plicit form, K\/2{z) — \fiTe~ z /\2z, which replaced in 
Eq. (|14|l and performing the resulting sum gives 

m\L) =a (r-r ) + ^ V |ilo g (l-e-^) 
+ OT [l ° g(1 -^ lL)]2 - 
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Taking m = 1 in Eq. l|16[l . we get the critical temperature, 

A plot of the relative critical temperature T c /Tq 
against the inverse thickness of the system from the 
above equation, is given in Fig. 2 for, in natural units, 
a = lm _1 ,/3 = 0.1m- 1 and T = 4000m- 1 (« 10°K). We 
see from the figure that the critical temperature slightly 
grows above the bulk transition temperature as the thick- 
ness of the system diminishes, reaching a maximum and 
afterwards starting to decrease until a zero value, cor- 
responding to a minimal allowed thickness for the sys- 
tem. This behaviour may be contrasted with the lin- 
ear decreasing of T c from the maximum value To with 
the inverse of the thickness of the system, that has been 
found for second-order transitions [0| UjJ ■ We also re- 
mark that in D = 3, for second-order transitions, one 
considers to = and that leads to the need of a pole- 
subtraction procedure for the mass In our case such 
a procedure is not necessary, as a first-order transition 
must occur for a non-zero value of the mass. This fact, 
together with the closed formula for the Bessel function 
for D = 3, allows us to obtain the exact expression 118|) 
for the critical temperature. 
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